455 


RENI 


ON THE APPLICATION OF THE CALCULUS OF 
QUATERNIONS TO THE THEORY OF THE MOON* 


Communicated 14 June 1847. 
[Proc. Roy. Irish Acad. vol. 1m (1847), pp. 507-20.] 


Sir William Rowan Hamilton made a communication respecting the application of the 
Calculus of Quaternions to the Theory of the Moon. 


I. At two Meetings of the Royal Irish Academy, in the month of July 1845, Sir William 
Rowan Hamilton had exhibited} and illustrated the following general equation of motion of 


a system of bodies, with masses m, m’, ..., and with vectors «,a’,..., and attracting each other 
according to Newton’s law: a ait 
— B i; =., 1 
di? ~ (a—a wet (a—a’)*} (1) 


He had, at that time, deduced from this equation the known laws of the centre of gravity, of 
areas, and of living force, for any such multiple system; and had shewn that the corresponding, 
but less general, equation of relative motion of a binary system, which (by changing œ — æ’ to a, 
and m+’ to M) becomes d?a M 


de ot a]( — a?) — a)’ (2) 


can be rigorously integrated by the processes of his new calculus of quaternions, so as to 
conduct, with facility, when the principles and plan have been caught, to the known laws of 
elliptic, parabolic, or hyperbolic motion of one of the two attracting bodies about the other. 
(See the Proceedings of 14 and 21 July 1845, Appendix to Volume ITI, pp. xxxvii, &c.)t 

At a subsequent Meeting of the Academy, in December 1845,§ Sir W. Hamilton had shewn 
that the general differential equation (1) might be put under this other form: 


d? a d2a 
di t dt? 


mm 
Vir a-e (3) 


and that it might, theoretically, be integrated by an adaptation of that ‘General Method in 
Dynamics’ which he had previously published in the Philosophical Transactions of the Royal 
Society of London, for the years 1834 and 1835;|| and which depended on a peculiar combina- 
tion of the principles of variations and partial differentials, already illustrated by him, in 
earlier years, for the case of mathematical optics, in the Transactions of this Academy. (See 
Proceedings of 8 December 1845, Appendix already cited, pp. lii, &c.)+ 

At the same Meeting of December 1845, Sir W. Hamilton assigned the two following rigorous 
differential equations for the internal motions of a system of three bodies, with masses m, m’, m”, 


0=22. m (0x a) +92- 


* [See also Elements, article 421.] + [See XLI.] 
t [See XLI.] ‘*§ [See XLII.] || [See XLIT, p. 450, footnote.] 
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and with vectors a, f + a, y+a—that is, for the motions of the two latter of these three bodies 
(regarded as points) about the former—as consequences of the general equation (1): 


ap m+m' sinea A aibs y 
wees EUA is Ion F 
dy m+m" ya m {2 i Ai -Eal 
d? -y4/(—Y*) J-0 -8 ie Py 

It was remarked, that by regarding m, m’, m”, as representing respectively the masses of 
the earth, moon, and sun, £ and y become the geocentric vectors of the two latter bodies; and 
that thus the laws of the disturbed motion of our satellite are contained in the two equations 
(4) and (5)—but especially in the first of those equations (the second serving chiefly to express 
the laws of the sun’s relative motion). 

The part of this equation (4), which is independent of the sun’s mass m”, is of the form (2), 
and contains the laws of the undisturbed elliptic motion of the moon; the remainder is the 
disturbing part of the equation, and contains the laws of the chief lunar perturbations. A com- 
mencement was made of the development of this disturbing part, according to ascending 
powers of the vector of the moon, and descending powers of the vector of the sun; and an 
approximate expression was myn obtained, which may be written thus: 


EA Wont a vo I> ame P+ 3YPY) (6) 
A-B} (= Cs ae La 

There was also given a geometrical CAM Ey of this result, corresponding to a certain 
decomposition of the sun’s disturbing force into two others, of which the greater is triple of 
the less, while the angle between them is bisected by the geocentric vector of the sun; and the 
lesser of these two component forces is in the direction of the moon’s geocentric vector 
prolonged, so that it is an ablatitious force, which was shewn to be one of nearly constant 
amount. 

Although the foregoing formulae may be found in the Appendices already cited,* to the 
Proceedings of the above-mentioned dates, yet it is hoped that, in consideration of the import- 
ance and difficulty of the subject, and the novelty of the processes employed, the Academy 
will not be displeased at having had this brief recapitulation laid before them, as preparatory 
to a sketch of some additional developments and applications of the same general view, which 
have since been made by the author. It may, for the same reason, be not improper here to 
state again, what was stated on former occasions, that all expressions involving vectors, 
a, a’, &c., such as are considered in this new sort of algebraical geometry, and enter into the 
foregoing equations, admit of being translated into others, which shall involve, instead of those 
vectors, three times as many rectangular coordinates, x, y, z, x’, y’,z’, &c., be means of relations 


of the forms a=ix+jy+ kz, a’ =ix' +jy' +kz', &c.; (7) 
where tj k are the three original and coordinate vector units of Sir William Hamilton’s theory 
of quaternions, and satisfy the fundamental equations 

P=opP=e=1; yY=k, jk=i, ki=j; jik, kj=-t, tkh=-37; (8) 
which were communicated to the Royal Irish Academy at the Meeting of 13 November 1843. 
(See the Proceedings of that date, and the author’s First Series of Researches respecting 


Quaternions, which Series has lately been printed in the Transactions of the Academy, Vol. XXI, 
Part 2.)+ 


(5) 


* [See XLI, XLII] t [See V, VIL] 
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II. It is evident, from inspection of the equations above recapitulated, that every trans- 
formation of the vector function, (a) =a-1(—a2)-4 (9) 


which represents, in direction and amount, the attraction exerted by one mass-unit, situated 
at the beginning of the vector æ, on another mass-unit situated at the end of that vector, must 
be important in the theory of the Moon; and generally in the investigation, by quaternions, of 
the mathematical consequences of the Newtonian Law of Attraction. The integration of the 
equation of motion (2) of a binary system was deduced, in the communication of July 1845,* 
from a transformation of that vector function, which may now be written thus: 
i 4 2d.a(—a?)t. 
AT | (10) 
where d is, as in former equations, the characteristic of differentiation. And the hodographic 
theory of the motion of a system of bodies, attracting each other according to the same 
Newtonian law, so far as it was symbolically stated to the Academy, at the meeting of 
14 December 1846, depends essentially on the same transformation. In fact, if we make 


da =r di, a= frat: (11) 


and if, by the use of notations explained in former communications, we employ the letters 
U and V as the characteristics of the operations of taking the versor and the vector of a 
quaternion, writing, therefore, 

U(a)=a(—a2)-#; V.er=—V.1a=}h(ar—72); (12) 


the equation (2) of the internal motion of a binary system becomes 


eau frat) 


V(r{rau) } 


where the denominator in the second member is constant, by the law of the equable description 
of areas. Hence, this second member, like the first, is an exact differential; and an immediate 
integration, introducing an arbitrary but constant vector e, coplanar with æ and 7, gives the 
law of the circular hodograph, under the symbolical form 


(13) 


vei alias) (14) 


Vrf dt 


the constant part of this expression (14) for the vector of the velocity, 7, being the vector of 
the centre of the hodograph, drawn from that one of the two bodies which is regarded as the 
centre of force; while the variable part of the same expression for 7 represents the variable 
radius of the same hodographic circle, or the vector of a point on its circumference, drawn 
from its own centre of figure as the origin. 


Multiplying this integral equation (14) by f dt, taking the vector part of the product, 
* [See XLI.] 


t [See Mathematical Papers, vol. 1, paper X, and Lectures, article 621 (for historical remarks on the 
Hodograph).] 
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dividing by M, and multiplying both members of the result into the constant denominator of 
the second member of (13) or of (14), we find, by the rules of the present calculus, 


2 
-(v.rfrar] 
Aea rdt+T. |r dt; (15) 


where S and T are the characteristics of the operations of taking respectively the scalar and 
tensor of a quaternion, so that, as applied to the present question, they give the results, 


T.[rdt=Ta= y(-a)=r; (16) 
and S.efr dt = 4(ex + ae) = er cos v; (17) 
where e=Te= ./(—e?) = const.; (18) 


while v is the angle (of true anomaly) which the variable vector a of the orbit makes with the 
fixed vector —e in the plane of that orbit; and r denotes the length of a, or what is usually 
called (and may still in this theory be named) the radius vector of the relative orbit. The first 
member of the equation (15) is a positive and constant number, representing the quotient 
which is obtained when the square of the double areal velocity in the relative orbit is divided 
by the sum of the two masses; if then we denote, as usual, this constant quotient (or semi- 
parameter) by p, and observe that the constant e is also numerical (expressing, as usual, the 
eccentricity of the orbit), we shall obtain again, by this process, as by that of July 1845,* 
the polar equation of the orbit, under the well-known form, 
p 

"=I Fecosv (19) 
This sketch of a process for employing the general transformation (10) in the theory of a binary 
system, may make it easier, than it would otherwise be, to understand how the following 
equation for the motion of a multiple system, 


(m+Am) dU ( far at) 


v (ar. farat) i 


(where m+Am, T+Ar7, are the mass and the vector of velocity of an attracting body, as 
m,T are those of an attracted one, which is analogous to, and includes, the equation (13) for 
the motion of a binary one, and which agrees with a formula communicated} to the Academy 
in December 1846), was obtained by the present author; and how it may hereafter be applied. 


dr == (20) 


III. The vector function ¢(«) in (9) may be called the TRACTOR corresponding to the vector 
of position a, or simply the tractor of æ; and another general transformation of this tractor, 
which is more intimately connected than the foregoing with the problem of perturbation, may 
be obtained by supposing the vector æ to receive any small but finite increment /; representing 
a new but shorter vector, which begins, or is conceived to be drawn, in any arbitrary direction, 
from the point of space where the vector « ends; and, by then developing, in conformity with 
the rules of quaternions, the new tractor 6(f+ a), (answering to the new vector +a, which is 


* [See XLI.] t [See XLIII, p. 457, second footnote.] 
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drawn from the beginning of æ to the end of £), according to the ascending powers of this added 
vector f. In this manner we find 


p(B +a)={-(B+a)}4(B+a)4 
={—a?(1 +2718) (1+ pa) {a(l+a1B)} 4 


= (1+ pa) (1+ a-18)-Fa(— a?) -4; (21) 

that is, pe + a) = B: n’ Pn, n’ (22) 
if we make, for abridgment, 

Pn, n= Mn, (pa) (xp) a — at- (23) 

A i Sha 3...(2n—1) 3.5... (2n'+1) (24) 


2.4...(2n) ~ 2.4...(2n’) ° 


The attraction ¢(#+a) which a mass-unit, situated at the beginning of the vector #+a, 
exerts on another mass-unit situated at the end of that vector, is thus decomposed into an 
infinite but convergent series of other forces, Øp w, of which the intensities are determined by 
the tensors, and of which the directions are determined by the versors, of the expressions included 
in the formula (23); er by the following expressions, which are derived from it by the rules of 
the calculus of quaternions: 


Tón w= ny (TE) (Lays (25) 


Ugn w= (Ufa) (Uaj+= (UÉ) "U -a). (26) 


Let a, b, be the lengths (or tensors) of the vectors «æ, 2, and let C be the angle between them, 
which angle we may conceive to express the amount of the positive rotation, in their common 
plane, from the direction of — æ to the direction of +£; then the positive or negative rotation 
in the same plane, from the same direction of —«, to the direction of the component force 
n,n» is expressed as follows: 


angle, from —«a to force ¢, w =(n—n')C; (27) 
n+n i 
and intensity of same component force = M, w £) a2, (28) 


The case n=0, n’=0, answers to the old tractor d(«), or to a force of which the intensity is 
represented by a~?, while its direction is the same as that of — æ. 


IV. Thus, if the vector æ be conceived to begin at a point B, and to end at the point C, while 
_the vector £ shall be conceived to begin at C, and to end at A; and if we conceive a unit-mass 
at B to attract two other masses, regarded as collected into points, and as situated respectively 
at C and at A; this attraction of B will disturb the relative motion of A about C, if A be supposed 
to be nearer than B is to C, by producing a series of groups of smaller and smaller forces, of which 
groups it may be sufficient here to consider the two following. 

The first and principal group consists of the two disturbing forces ¢, o and ġo, ,, and of these 
the first is purely ablatitious, or is directed along the prolongation of the side of the triangle 
ABC, which is drawn from C to A, and it has its intensity denoted by the expression }ba-%, 
since we have for this force, and for its tensor and versor, the expressions 


10=2h(—2) 4; Td, 9=3ba*; Ud, .=Uf. (29) 
The second disturbing force, of this first group, has for expression 
$o,1= apa — a?) = apaa; (30) 
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its intensity is exactly triple of that of the former force, being represented by 3ba-%; and its 
direction is the same as that of a straight line drawn from C to A’, if A’ be a point such that 
the line AA’ is perpendicularly bisected by the line BC (prolonged through C if necessary). 
These two principal disturbing forces evidently correspond to those which were considered 
for the case of our own satellite in a communication above alluded to; the second force being 
the one which was described in that former communication as being directed to what was 
there called the ‘fictitious moon,’ and was conceived to be as far from the sun in the heavens 
on one side, as the actual moon is on the other side, but in the same great circle. 

If we now extend that mode of speaking so far as to conceive a similar reflexion of the sun 
with respect to the moon, and to call the point in the heavens so found the ‘fictitious sun,’ the 
moon being thus imagined to be seen midway among the stars between the actual and the 
fictitious sun: and if we farther imagine a ‘second fictitious sun,’ so placed that the actual 
sun shall appear to be midway between this and the first fictitious sun; we shall then be able 
to describe in words the directions of the three disturbing forces of the second group, and to 
say that they tend respectively, for the case of our own satellite, to these three (real or 
fictitious) suns. For these three forces will have, for their respective expressions, the three 
corresponding terms of the development of the tractor (22), namely, the following: 


$o,o= 8haP(—a%)4; $1 1=3f%a(-a%)4; $4 .=1afapfa-\—a%)4; (31) 
of which the intensities are respectively 
gbat; #b7a-*; 42b7a-4; (32) 


so that they are exactly proportional to the three whole numbers, 1, 2, 5; while they are directed, 
respectively, to the first fictitious sun, the actual sun, and the second fictitious sun. The 
disturbing force of a superior planet, exerted on an inferior one, may be developed or decom- 
posed into a series of groups of lesser disturbing forces, the intensities of the several forces in 
each group being constantly proportional to whole numbers, in an exactly similar way; nor 
does the application of the principle and method of development thus employed terminate 
here. In the applications to the lunar theory, a and b, in the recent expressions, are to be 
regarded as denoting the variable distances of the sun and moon from the earth; and the 
expressions for the forces are to be multiplied by the mass of the sun. Nothing depends, so far, 
on any smallness of eccentricities or inclinations. 


V. The lunar theory is, very approximately, contained in the differential equation (4), 
provided that we regard y as the elliptic vector of the sun, drawn from the common centre of 
gravity of the earth and moon; and the laws of the sun’s relative elliptic motion, with respect 
to that centre of gravity, are then contained in the following differential equation, which takes 


the place of the equation (5): E an ec Hee 


ae a cela 
Indeed, when we come to consider the small disturbing forces which belong to the second 


group, and which depend on the inverse fourth power of the sun’s distance, the corresponding 
terms of the development of the first member of the formula (6) are then, for greater accuracy, 


(33) 


to be multiplied by the fraction ae , which expresses the ratio of the difference to the sum 


of the masses of the earth and moon. But if we neglect, for the present, those small disturbing 
terms, we may regard that formula (6) as accurate, without as yet neglecting anything on 
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account of smallness of eccentricities or of inclinations; and even without assuming any 
knowledge of the smallness of the moon’s mass, as compared with the mass of the earth; y still 
denoting, as just stated, the elliptic vector of the sun. And thus, if the moon’s geocentric 
vector J be changed to the sum 2 +82, where the term df is supposed to depend on the dis- 
turbing force, and to give a product which may be neglected when it is multiplied by or into 
the expression for that force, we shall have the following approximate differential equation, 
by developing the disturbed or altered tractor ¢(f+6f), and confining ourselves to the first 
power of df: 


d?ôf d? m” 
ee le aie ra= s. pÈ 5 (f+ 3A-68 6) + 5s 2(— y2)8 (B+3y—"By). (34) 


The disturbance ôf of the moon’s geocentric vector is thus exhibited as giving rise to an 
alteration 6¢(f) in the corresponding tractor $(f), which alteration is analogous to a disturbing 
force, and occasions the presence of the first of the two parts of the second member of the 
equation (34): which equation will be found to contain a considerable portion of the theory 
of the moon. 


VI. The author will only mention here two very simple applications, which he has made of 
this equation (34), one to the Lunar Variation, and the other to the Regression of the Node. 
Treating here the sun’s relative orbit as exactly circular, and the moon’s as approximately 
such, neglecting the inclination, taking for units of their kinds the sum of the masses of the 
earth and moon, and the moon’s mean distance and mean angular velocity, and employing, as 
usual, the letter m to denote (not now the earth’s mass, but) the ratio of the sun’s mean angular 
motion to the corresponding motion of the moon, the differential equation (34) becomes: 


d?ô, m? 
OR = HOB + 30968) +" (B+ 37y); (35) 
in which the laws of the circular revolutions of the vectors 2 and y give 
d? d ; 
p; Y=- my. (36) 


Assuming, from some general indications of this theory, an expression for the perturbation 
of the moon’s vector, which shall be of the form 


dB =m (AP + By py + Ch y"Byp), (37) 

and neglecting all powers of m above the square, we find 
TPP = — m? Af + Bytpy + 3°0B-1y-16yf); (38) 
B68. B= m(AB + Cy Py + BBy"Byp); (39) 
so that the three numerical coefficients, A, B, C, must satisfy the three following equations of 
condition: —~A=2A+}; giving A=-}; (40) 
and : —~B=}(B+3C)+3; —90=}(C+3B); (41) 
giving ` B=-23; C=+,%. (42) 


Thus, if we neglect eccentricities and inclinations, and confine ourselves to the first power 
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of the disturbing force, or to the second power of m, the perturbation of the moon’s vector, 
produced by the sun’s attraction, is composed of the three following terms: 


p= -E p- apy 4 2 pay-apyp. (43) 


The first of these three terms expresses that the sun’s ablatitious force, by partially counter- 
acting the earth’s attractive force on the moon, allows our satellite to revolve in a somewhat 
smaller orbit than would otherwise be consistent with the observed periodic time: the ratio of 

2 
the diminished to the undiminished radius of the orbit being that of 1 -7 to 1. The second 
term expresses a displacement of the moon, through perturbation, from its diminished circular 
orbit, of which displacement the constant magnitude or length bears to the radius of the 


ne y . p 19m? í 5 EAN AS 
undiminished orbit the ratio of ien to unity; while the direction of this displacement is always 


from that fictitious moon, to which it has been seen that one of the two principal components 
of the sun’s disturbing force is directed: an opposition of sign which may at first surprise, but 
which is exactly analogous to the contraction of the orbit produced by the ablatitious force 
(when the periodic time is given), and is to be explained upon similar principles. Finally, the 
third term of the formula (43) for ôf, expresses that with the two foregoing displacements 
a third is to be combined, which is, like them, of constant amount, being equal to ths of 


2 
the second displacement, or bearing to the radius of the moon’s orbit the ratio of pine to unity; 


but being always directed to what, by an extension of a recently employed phraseology, might 
be called the second fictitious moon, being so placed that the actual moon is midway in the 
heavens between this fictitious moon and the one which was before considered. These two 
latter terms of (43) contain the chief laws of the Lunar Variation; and are easily shown to give 
the known terms in the expressions of the moon’s parallax and longitude, 


Per ies. | (44; 45) 


8 =m* cos 2() — ey apati 


It may assist some readers to observe here, that when the inclination of the orbit is neglected, 
the longitudes of the first and second fictitious moons are, respectively, 


2©-), and 3)-—20; (46) 
while those of the first and second fictitious suns, mentioned in a former section of this abstract, 
are, under the same condition, ay 2G! andi ae lan: (47) 


VII. The law and quantity of the regression of the Moon’s Node may also be calculated on 
principles of the kind above stated, but we must content ourselves with writing here the 
formula for the angular velocity of a planet’s node generally, considered as cae on the 


da 
variable vector of position a, the vector of velocity ai? and the vector of acceleration — de’ * and also 


on a vector unit A, supposed to be directed towards the north nh of a fixed ecliptic. The 
formula thus referred to is the following: 


S.à. S.d? daq 


VN ai (48) 


dQ= 
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where S and V are, as before, the characteristics of the operations of taking the scalar and 
vector of a quaternion. The author proposes to give a fuller account of his investigations on 
this class of dynamical questions, when the Third Series of his Researches respecting Quater- 
nions shall come to be printed in the Transactions of the Academy: the Second Series being 
devoted to subjects more purely geometrical; as the First Series (already printed) relates 
chiefly to others which are of a more algebraical character.* 


* [After VII, Hamilton published no other work on Quaternions in the Transactions. See Elements, 
Book mm, Chapter m1. ] 
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